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We examine the effect of cluster size on the interaction of Ar 5 5-Ar2057 with intense extreme ultra- 
violet (XUV) pulses, using a model we developed earlier that includes ionization via collisional exci- 
tation as an intermediate step. We find that the dynamics of these irradiated clusters is dominated 
by collisions. Larger clusters are more highly collisional, produce higher charge states, and do so 
more rapidly than smaller clusters. Higher charge states produced via collisions are found to reduce 
the overall photon absorption, since charge states of Ar 2+ and higher are no longer photo- accessible. 
We call this mechanism collisionally reduced photoabsorption, and it decreases the effective cluster 
photoabsorption cross-section by more than 30% for Arss and 45% Ar2057. An investigation of the 
shell structure soon after the laser interaction shows an almost uniformly charged core with a mod- 
estly charged outer shell which evolves to a highly charged outer shell through collisions. This leads 
to the explosion of the outer positive shell and a slow expansion of the core, as was observed in mixed 
clusters at shorter wavelength [ ]. The time evolution of the electron kinetic energy distribution 
begins as a (mostly) Maxwellian distribution. Larger clusters initially have higher temperature, but 
are overtaken by smaller temperature after the laser pulse. The electron velocity distribution of 
large clusters quickly become isotropic while smaller clusters retain the inherent anisotropy created 
by photoionization. Lastly, the total electron kinetic energy distribution is integrated over the spa- 
tial profile of the laser and the log-normal distribution of cluster size for comparison with a recent 
experiment [2], and good agreement is found. 



I. INTRODUCTION 

The interaction of atomic clusters with intense ultra- 
fast laser pulses has been investigated over a range of 
wavelengths in recent years. Most of this work has been 
in the infrared (IR) regime [ ], but as new very intense 
shorter wavelength sources have come online, this inter- 
action has been investigated up to the soft-Xray [1, 4-10]. 

Cluster interaction with intense laser pulses is highly 
wavelength dependent. In the IR, tunnel ionization and 
electron heating processes dominate and the subsequent 
plasma dynamics are driven by the laser field itself [3, 11]. 
This is not true at shorter wavelengths. Experiments on 
clusters in the vacuum ultraviolet (VUV) regime near 100 
nm [12] observed unexpectedly high ionic charge states 
and a much larger energy absorption than predicted by 
existing models [8] . This sparked a concerted theoretical 
effort over the last eight years. Santra et al. proposed 
an enhancement to inverse bremsstrahlung heating (IBH) 
due mainly to using self-consistent potentials for ions [13] 
and later plasma screening effects [14]. Jungreuthmayer 
et al. proposed a many-body dielectric recombination 
scheme for the strongly-coupled plasma electrons in the 
cluster whereby the electrons could then be driven by 
the laser field to reionize [15]. Siedschlag et al. proposed 
allowing photoionization to occur above the classical po- 
tential barrier of the neighboring ion [10]. Ziaja et al. 
incorporated many of these and found they all played 
various roles depending on the intensity [ ]. It is clear 
from these works that clusters irradiated by short laser 
pulses with wavelengths in regime of single photon ion- 
ization represent a new and theoretically challenging area 
of physics. 



Recent experiments around 40 nm have begun to probe 
the regime in which the photoelectron has a significant 
amount of kinetic energy and ions are photoaccessible 
[2, 18]. Contrary to the shorter wavelength regimes, the 
laser field coupled very weakly to the cluster electrons, 
producing a photoelectron spectrum with almost no sig- 
nal above the atomic photoelectron energy in gas. The 
lower energy region showed that the cluster was charg- 
ing and cooling subsequent photoelectrons down to lower 
energy. The clusters also produced high charge states 
well above what was detected in gas. This interaction 
was also found to depend on cluster size [2, 18]. Fur- 
ther, it was found that even small clusters emitted ions 
with very little kinetic energy (less than 30 eV) [2]. An 
attempt to address the charge transfer and explore the 
explosion dynamics was done using heterogeneous clus- 
ters [1, 19]. Many aspects of these experiments are not 
well understood even in gases [20]. 

The XUV regime offers unique opportunities in the 
study of the dynamics of finite systems such as clusters. 
Single photon ionization from the ground state of rare gas 
atoms is accessible, though inner shell ionization is not. 
In addition, unlike in the IR and VUV, the IBH of freed 
electrons is negligible, due to the very low quiver energy, 
decoupling cluster plasma dynamics from the laser field 
[ ]. However, electrons ionized from atoms have a sig- 
nificant amount of kinetic energy, enough to collisionally 
ionize neutral Argon. To date, there have only few the- 
oretical works in this regime, however. Bodstedt et al. 
proposed a multistep model to explain the photoelectron 
spectrum of their experiment on small Argon clusters ex- 
posed to intense 32 nm laser pulses [2] . Single photon ion- 
ization events were determined based on a Monte-Carlo 
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model, and after each ionization step an electron is emit- 
ted from the cluster with an energy that depended upon 
the building positive space charge of the cluster. Ziaja et 
al. also examined the photoelectron emission spectrum, 
but used a kinetic Boltzmann equation technique that 
was then compared to molecular dynamics simulations 
results [ ]. They found good agreement with the mea- 
sured electron emission spectrum and showed that IBH 
was insignificant at 32 nm. Arbiter et al. used a molec- 
ular dynamics and Monte Carlo technique to look at the 
electron emission spectrum at different intensities finding 
that at high intensity there are a significant number of 
thermal electrons emitted from the cluster [22]. 

Most recently, we showed the importance of two-step 
collisional ionization via an intermediate excited state 
[23]. This model allowed atom and ions to first be ex- 
cited by an impact electron, which requires much less 
energy than ionization, and then ionized from the ex- 
cited state. We called this model augmented collisional 
ionization (ACI) and were able to reproduce the high- 
est charge state seen in Ref. [ ]. This paper will expand 
upon our previous work to give a more detailed analy- 
sis of Argon cluster interaction with intense 32 nm laser 
pulses, including examining the cluster explosion dynam- 
ics and how these depend upon cluster size. Understand- 
ing cluster dynamics in the XUV regime, while inter- 
esting in its own right, may also inform planned single 
shot large molecule imaging experiments at even shorter 
wavelengths [24, 25]. 

The paper is organized as follows. In section II we de- 
scribe the model and including details of how photoion- 
ization, collisional ionization and augmented collisional 
ionization (ACI) are implemented. Section III presents 
the results for Argon clusters of size 55, 147, 561 and 2057 
atoms, with Section III A concentrating on ions and Sec- 
tion IIIB on electrons. In Section IIIB3 we obtain the 
calculated photoelectron spectrum corresponding to the 
experimental results of Ref. [ ], and find good agreement. 
Finally, Section IV gives a summary of all the results. 
Atomic units (H = e = m e = 1) are used throughout 
unless otherwise specified. 



II. THEORY 

In this work we use hybrid approach to simulate the 
time evolution of Argon clusters in an XUV radiation 
field, where classical molecular dynamics is used for ion 
and electron motion and a Monte-Carlo scheme is used 
to determine ionization events. Ionization rates are de- 
termined by the quantum mechanical transition cross- 
section for the various processes. The simulation begins 
with a neutral cluster of atoms. Ionization events re- 
sult in electrons being created within the code, and in an 
increment in the charge state of the the parent neutral 
or ion. Collisional excitation events are also permitted, 
whereby a neutral or ion may become excited; in this case 
no new electron is added to the simulation but the par- 



ent neutral or ion is set to be in an excited state, which 
determines its future ionization potential and ionization 
cross-section. 

In order to mitigate numerical heating, the classical 
motion of particles was calculated assuming the particles 
were Gaussian distributions, instead of point particles, 
via the following smoothed potential 

{Q/r r > a 

QSexp(-I(^) 2 ) r<* W 

where r is the radial distance, B is the maximum poten- 
tial depth for a singly-charged ion, Q is the charge, and 
cr is an effective smoothing radius given by 

-|exp(i). (2) 

Note that (j> is continuously diffcrentiable at the transi- 
tion point r — a. To ensure energy conservation, the 
potential of an electron and a singly-charged ion are set 
to be equal in magnitude, and the potential depth of the 
multiply charged ions are integer multiples of that of the 
singly-charged ion's potential depth. Newly ionized elec- 
trons are created at the same location as their parent ion 
to avoid dipole heating. 

Ionization in the cluster environment is modeled as an 
isolated system within a background cluster potential. 
We take the cluster potential as constant over the outer- 
most electron's wavefunction. Thus the cluster potential 
at the location of an ion or neutral is taken as the thresh- 
old for ionization for that ion or neutral. All ionization 
processes are then calculated with respect to this thresh- 
old. This allows for the use of atomic photoionization 
and collisional excitation/ionization cross-sections with 
only a small error due to the approximation [23]. 

In the XUV, direct photoionization is the only way 
for the radiation field to deposit energy into the cluster. 
The single photon ionization probability is determined at 
each time-iteration for every neutral/ion, depending on 
its charge state and the photon flux. The photon flux 
is determined via the intensity, which is modulated both 
by the time profile of the pulse and by photon absorp- 
tion. After each photoionization event the intensity of 
the pulse is decreased by one photon. This more accu- 
rately models low flucncc pulses. The cross-section for 
the photoionization of neutral Argon was obtained from 
Ref. [26]. Those of ionic Argon were obtained using Los 
Alamos Atomic Physics Codes [27]. Multiphoton ioniza- 
tion is negligible at this photon fluence [2]. 

Two channels are available for collisional ionization: a 
single-step transition from the ground state and a two- 
step transition through an intermediate excited state 
called augmented collisional ionization (ACI) [ ]. The 
cross-section of all energetically accessible states (excita- 
tion and ionization) are combined to give the total cross- 
section for the occurrence of an event (whether excita- 
tion or ionization). If the impacting electron's trajectory 
is within this cross-section an event will take place. A 



3 



Monte Carlo scheme is used to determine which type of 
event takes place based on the relative weights of the 
cross-section for each transition. Single-step collisional 
ionization cross-sections were calculated using the semi- 
empirical Lotz formula and coefficients [28]. The Born 
plane-wave approximation [29] was used for both excita- 
tion and ionization from an excited state. 

The excited states considered consist of a subset of 
all possible excited states. Only single electron excited 
states were used and of those, only the lowest eight states 
with I < 3 were implemented. This subset is the most im- 
portant as it contains the lowest energy states. Including 
more single electron excited states adds to the total col- 
lisional cross-section, although states near the threshold 
require almost as much energy as ionization and are thus 
almost as infrequent. The energies of the excited states 
and the cross-sections were obtained using the Hartree- 
Fock implementation of the Cowan code [29]. 

III. RESULTS 

In this section we illuminate the details of the laser 
cluster interaction and determine the effects of cluster 
size on key measures of the cluster dynamics. These 
include: ion charge state evolution and distribution, 
excited state evolution, the relative importance of the 
relevant ionization channels, and electron energies and 
charge transfer. Generally, the interaction dynamics pro- 
ceeds as follows. After the first few direct photoioniza- 
tions, the cluster builds up a positive space charge such 
that at some point subsequent photoelectrons are pre- 
vented from escaping the cluster; this is described Ref. [2] 
via their multistep ionization model. Our model further 
allows the cluster-bound photoelectrons to either release 
further electrons via collisional ionization or cause colli- 
sional excitation of neutrals or ions in the cluster during 
and after the laser pulse. 

The parameters of the radiation field for each data set 
were A = 32.8 nm, I = 5 x 10 13 W/cm 2 and a full-width 
at half-maximum of 25 fs. Closed-shell icosahedral Ar- 
gon clusters of 55 (2 shells), 147 (3 shells), 561 (5 shells) 
and 2057 (8 shells) were used. These clusters were re- 
laxed according to a Lennard- Jones potential for neutral 
Argon. The results presented are for clusters assumed to 
be in the laser focus, and each plot is an ensemble average 
over many simulations, the number of which was chosen 
to ensure that more than 4 x 10 atoms are included in 
the average. 

A. Ions 

The charge states of the ions are a signature of the dy- 
namics of intense laser cluster interaction. Charge states 
higher than what is accessible through photoionization is 
an indication that there are cluster-driven processes that 
are not present in intense laser interaction with the gas 



phase of the material. This occurs in clusters when the 
photon fluence is high enough to have at least a mod- 
est amount of photoionization within a single cluster, so 
that collisional processes become relevant. The parame- 
ter range of this work is such that the interaction is well 
above the collisional threshold. 



1. Charge state 

In Argon gas targets, only the Ar 2+ is accessible via 
direct photoionization. However experiments with clus- 
ters have detected higher charge states indicating the 
importance of cluster-driven processes such as collisions 
[ ]. Figure 1 plots our calculated charge state distri- 
bution for Argon clusters of sizes (a) 2057, (b) 561, (c) 
147 and (d) 55, for the laser parameters listed above. 
Though photoionization stops after the laser pulse has 
passed, collisional ionization and excitation can still oc- 
cur. Thus, these snapshots were taken after these pro- 
cesses also stopped, at one picosecond after the pulse. 
The bare charge states are reported. 

We see that the average charge state increases with 
cluster size. This is despite the fact that larger clus- 
ters have fewer photons per atom compared with smaller 
clusters. However, larger clusters develop a larger space 
charge, and thus there are a greater number of photoelec- 
trons that become bound to the cluster. These can then 
precipitate a collisional ionization or excitation event. 
Further, in a larger cluster there are more available colli- 
sional ionization or excitation targets. As we will see later 
in section III A 4, collisional processes dominate over di- 
rect photoionization in charge state creation above Ar 1+ . 

Figure 2 gives the yields per atom of Ar 3+ (blue dots) 
and Ar 4+ (red solid line) as a function of cluster size. Ad- 
ditional cluster sizes were considered for this graph. The 
yields of both the Ar 3+ and Ar 4+ increase with cluster 
size, and moreover, the Ar 4+ yield increases more rapidly 
than the Ar 3+ . This rapid increase of the Ar + is thus 
not explained simply as an increase in Ar 3+ , but indi- 
cates that collisional processes become more important 
for larger clusters. This is because it is only through col- 
lisional processes that charge states above Ar 2+ can be 
created. 



2. Charge state evolution 

Our model also allows for the detailed tracking of 
each ion species during the interaction. The normal- 
ized population of each charge state as a function of 
time is shown in Fig. 3 for different sized clusters. Each 
plot shows the fraction of atoms that are neutral (blue 
dashed-dotted lines), Ar 1+ (green medium-dashed lines), 
Ar 2+ (red sparsely-dashed lines), Ar 3+ (cyan long and 
short-dashed lines), Ar 4+ (magenta solid lines) and Ar 5+ 
(yellow short-dashed line). At the bottom of the figure, 
we include a plot of the temporal profile of the Gaussian 




FIG. 1. (color online) The charge species yield per atom for: 
(a) Ar2057, (b) Ars6i, (c) Ari.47 and (d) Arss clusters interact- 
ing with a 25 fs, 32 nm laser pulse of 5xl0 13 W/cm 2 . These 
snapshots were taken one picosecond after the pulse, after 
collisional events were no longer occurring. 



laser pulse. Note that a logarithmic scale is used for the 
time axis. 

All clusters begin neutral and are initially ionized pri- 
marily by photoionization. The behavior of the low 
charge species early on is qualitatively similar for all clus- 
ter sizes. The population of neutrals decreases rapidly 
and is surpassed by the Ar 1+ followed by the Ar 2+ . The 
Ar 2+ remains the most abundant charge state for the two 
smaller clusters. The larger clusters continue to have ion- 
ization well past 100 fs, and at around 300 fs the Ar 3+ 
becomes (and remains) the most abundant. After 1 ps, 
almost no further collisional ionization occurs. 
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Cluster size 

FIG. 2. (color online) The yield per atom of the Ar 3+ is shown 
by the (blue) dotted line and the Ar 4+ by the solid (red) line 
as a function of cluster size. The shaded (gray) region is only 
an aid to view how the difference between the lines diminishes 
with cluster size. 
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FIG. 3. (color online) The charge species yield per atom as a 
function of time for: (a) Ar2057, (b) Arsei, (c) Ari.47 and (d) 
Ar 55 clusters interacting with a 25 fs laser pulse of 5xl0 13 
W/cm 2 . The bottom plot is the temporal profile of the 25 fs 
laser pulse. 
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The abundance of the higher charge states is quite 
small during the actual laser pulse. Contrary to longer 
wavelength regimes which are dominated by IBH, here 
photoionization is the only direct contribution of the 
laser pulse to charge state creation. The absorption of 
the XUV photons starts the ionization process, but the 
higher charge states appear in significant numbers only 
after the cluster has absorbed a sufficient amount of en- 
ergy (as is further demonstrated in Section IIIB1). At 
the intensity we consider, this begins only near the tail 
end of the laser pulse and continues for several hundred 
femtoseconds. 

The crossings of the curves in Fig. 3 occur earlier for 
larger clusters. Given that there are less photons per 
atom for the larger clusters, this cannot be due to in- 
creased photoionization, and thus indicates that larger 
clusters are more efficient at dispersing the energy from 
the laser through increased collisions. 

The location of the low charge state crossings all occur 
about 7 fs earlier for the Ar 2 o57 compared with the Ar 55 . 
This is almost 30% of the full- width at half-maximum of 
the laser pulse, and results in a significant decrease in 
how much energy the larger cluster can absorb from the 
laser. This is because more ions become transparent to 
the laser earlier due to collisions in the larger cluster, and 
thus do not absorb as many photons. If photoionization 
were the dominant process, the low charge state crossings 
for larger clusters would occur later, not earlier, than for 
the smaller clusters, due to the larger number of photons 
per atom for the smaller clusters. 

The fact that the neutrals decrease in population faster 
in larger clusters due to increased collisional ionization 
processes has other consequences. Neutrals have the 
largest photoionization cross-section. As their popula- 
tion decreases the overall photoionization cross-section of 
the cluster as a whole will also decrease, further decreas- 
ing the importance of photoionization in the larger clus- 
ters. Further, photoionization increases the overall tem- 
perature of cluster-bound electrons, while it is decreased 
by collisional ionization. Thus, the electron temperature 
would be smaller when collisional ionization of neutrals 
is included because cluster-bound photoelectrons ionized 
from the neutrals contribute more energy than those pho- 
toionized from the Ar 1+ . Thus a decrease in number of 
photoelectrons released from the neutrals would result 
in a decrease of the electron temperature. Thus if these 
effects are neglected, this could overestimate energy ab- 
sorption of the cluster. 

We quantify the reduction of total cross-section and 
electron temperature due to collisions by performing our 
simulations with and without collisional processes in- 
cluded. Without collisional processes, we find Ar2057 ab- 
sorbs 1840 photons on average, 1401 by neutrals and 439 
by Ar 1+ . With collisional processes included, Ar 2 o57 ab- 
sorbs 1151 photons, 618 by neutrals and 533 by Ar 1+ . 
Thus the total number of photons absorbed is decreased 
by 37%, resulting in a 46% reduction in the amount of 
energy the electrons obtain from the laser. We call this 



collisionally reduced photoabsorbtion (CRP). 

The effect is smaller in the less collisional smaller clus- 
ters. Without any collisional processes Ar^ absorbs 50 
photons on average, 38 by neutrals and 12 by Ar 1+ . With 
all collisional processes included Ai^ absorbs 38 photons, 
23 by neutrals and 15 by Ar 1+ . Thus the total number of 
photons absorbed is decreased by 24% in Ar 55 clusters. 
This still results in a 31% reduction in the amount of en- 
ergy the electrons obtain from the laser. While smaller 
than for the Ar 2 Q57 clusters it is not a negligible effect. 



3. Excited states evolution 

The efficacy of collisional ionization is driven largely 
by the access to excited intermediate states. Figure 4 
gives the excited state yield per atom of all neutrals and 
ions corresponding to the simulations of Fig. 3. Also 
shown is the total number of excited ions as the solid 
(black) line. The curves are very similar in shape to 
Fig. 3. Close to 20% of all neutrals and ions for each 
charge species is excited at any given time. The pro- 
portion of excited species, for the most part, increases 
with the charge state. This is because the energy gap 
from an excited state to the threshold is larger for higher 
charge states. For lower charge states this energy gap 
is smaller, meaning that these excited species are com- 
paratively shorter lived. The lower energy electrons are 
more abundant than higher energy ones, and thus ex- 
cited states with a smaller gaps to threshold are more fre- 
quently ionized. The proportion of excited ions is higher 
at later times because collisions are much less frequent. 

The evolution of the ionic and excited state popula- 
tions show that the laser cluster interaction in the XUV is 
a process begun by photoionization but where collisional 
and photoionization processes occur in tandem, signif- 
icantly affecting each other as quantified in Section 3. 
The large amount of Ar 3+ in Fig. 3 and excited neutrals 
and ions shown in Fig. 4, where both snapshots are dur- 
ing the laser pulse, indicate that a significant amount 
of collisional ionization is occurring as these species can 
only be created via collisional processes. Further, the 
collisional processes modify the photoabsorbtion rate of 
the cluster. More excited neutrals (as found in the larger 
clusters in particular) indicate that more neutrals could 
eventually be collisionally ionized via augmented colli- 
sional ionization (ACI), which is the dominant cause of 
CRP as will be shown in the next section. 

Once sufficient energy has been deposited into the sys- 
tem by the laser, collisional processes (most notably ACI) 
disperse the energy throughout the cluster. The larger 
the cluster the more rapidly the cluster becomes dom- 
inated by collisional processes and thus reaches higher 
charge states. 
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FIG. 4. (color online) Number of excited neutrals and ions, 
normalized to the total number of atoms, as a function of 
time for: (a) Ar2057, (b) Ars6i, (c) Ari.47 and (d) Ars5 clusters 
interacting with a 25 fs laser pulse of 5xl0 13 W/cm 2 . The 
solid (black) line is the total number of excited neutrals and 
ions. 



4- Mechanisms of ionization 

In order to examine the influence of each ionization 
mechanism directly, we plot for each charge state in Fig. 5 
the percentage of that charge state population that was 
ionized by the various available ionization mechanisms. 
The diagonal-lined boxes (magenta) give the percentage 
ionized by photoionization. The diamond-filled boxes 
(red) give the percentage ionized by standard one-step 
collisional ionization. The filled boxes (blue) give the 
percentage ionized via ACL 

For 32 nm radiation, only the neutral and Ar 1+ ions 
can be photoionized. For all cluster sizes, photoioniza- 
tion of neutrals is proportionally larger than for the Ar 1+ , 
though it decreases for larger clusters, from 44% for Ar 5 5 
to 30% to Ar2057- The photoionization of the Ar 1+ de- 
creases from 36% for Ar 55 to 28% for Ar 2 o57- 

The photoionization cross-section of Ar 1+ is close to 
half that of the neutral [2] , though this ratio is not born 
out in Figure 5. This is further evidence of how col- 
lisional processes affect the dynamics. We find that a 
significant number of neutrals are being ionized by colli- 
sional processes, predominantly ACI, before photoioniza- 
tion can occur. There is thus a change due to CRP in the 
expected photoionization yield for the cluster compared 
with what is expected for the same number of atoms in 




FIG. 5. (color online) Percentage of events that led to the 
ionization of the indicated charge state, for all the relevant 
mechanisms, for (a) Ar2057, (b) Arsei, (c) Ari.47 and (d) Arss. 
The ionization mechanisms are: photoionization , single-step 
impact ionization and augmented collisional ionization shown 
by (magenta) diagonal lines, (red) diamonds and solid (blue) 
fills respectively. 

gas for the same photon fluence. 

The proportion of neutrals ionized by single step col- 
lisional ionization is very similar for all cluster sizes at 
around 17%. While it does in fact occur for higher charge 
states, it is below the 1% range (thus not visible in the 
graph) demonstrating that ACI dominates. The propor- 
tion of neutrals and Ar 1+ ionized by ACI is larger for 
larger clusters. This accounts for the relative drop in the 
photoionization since the single-step collisional ionization 
remains roughly constant. This is due to the system be- 
coming more collisional at earlier times for larger clusters 
increasing the effect of CRP. 

The roles of photoionization and collisional processes 
are thus very clear. In this regime of intensity and wave- 
length, the system is initially driven by photoionization 
but its evolution is shaped by collisional processes. The 
higher charge states appear in large numbers well after 
the pulse is over, created by collisional processes which 
become more important for larger clusters. 



5. Charged shell structure 

We now consider the spatial distribution of the charge 
states within the clusters. The initial icosahedral struc- 
ture has closed shells which remain largely intact due to 
the force on the ions being primarily from their mutual 
repulsion [ ] . Thus to understand the charge state distri- 
bution, we consider the net charge per atom as a function 
of shell index. The net charge within a shell is the to- 
tal charge of all particles in that shell including ions and 
electrons classically bound to the ions; this is then di- 
vided by the number of atoms in the shell to obtain the 
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FIG. 6. (color online) Net charge per atom versus shell index 
for (a) Ar 2 o57, (b) Ar S6 i, (c) Ari 47 and (d) Ar 55 as determined 
at the end of the simulations. The solid (blue) line is the 
average of the net charge per atom over all the shells. 



net charge per atom. This is plotted versus shell index 
in Fig. 6 for Ar2057, Ars6i, Ari47 and Ar^ clusters at 
the end of the simulation, i.e., after ionization processes 
have ceased. 

The average over all shells of the net charge per atom 
is shown in as the horizontal blue solid line, and it gives 
a measure of how the cluster as a whole is charged. It 
decreases with the increasing cluster size, indicating that 
while larger clusters access higher charge states the elec- 
trons remain bound to the cluster. Consistent with previ- 
ous work [1], the net charge per atom of the outer shells is 
found for all cluster sizes to be much higher than that of 
inner shells. Moreover, only the outermost shell is above 
the average. Therefore in all cases most of the charge 
resides on the outermost layer(s) of the cluster. Taking 
the quasi-neutral core to consist of those shells which are 
below the mean, we find that the core size generally in- 
creases with increasing cluster size, leaving the outermost 
shell to explode fastest. 

We can fit the net charge per atom versus shell number 
by an exponential function, 

f(a) = a(exp(6(s - 2)) - 1), (3) 



Cluster size 


a 


b 


2057 


0.041 


0.394 


561 


0.062 


0.562 


147 


0.157 


0.646 


55 


0.230 


0.789 



TABLE I. Fitting parameters for the net charge per atom 
versus shell index, according to the function given by Eq. 3. 



where s is the shell number, and a and b are fit parame- 
ters; b quantifies the disparity in mean charge state across 
the shells. We plot these fits as black lines in Fig. 6, and 
list the fitting parameters in Table I. The value of b in- 
creases with decreasing cluster size, showing that smaller 
clusters have a larger disparity in the net charge per atom 
between shells, particularly between the outermost and 
next outermost shells. Larger clusters have a smaller 
disparity, since the larger space charge retains more elec- 
trons overall. 

In Fig. 7 we plot a snapshot of the net charge per atom 
versus shell index shortly after the laser pulse (at 60 fs). 
The result is quite different from Fig. 6, which gives a 
snapshot long after the laser pulse when all collisional 
ionization processes have ceased. In Fig. 7 the net charge 
per atom is more evenly distributed over the shells. This 
implies that the shell structure is a consequence of charge 
migration through collisional ionization which is effec- 
tuated by the electrons. For larger clusters, however, 
charged outer shells are already emerging, whereas for 
the Arss and Ari47 clusters they are still quite homoge- 
neous. This, along with Fig. 3, suggests that the smaller- 
clusters take longer to become collisional. The forma- 
tion of charged outer shells early on in the larger clusters 
is another indication of a significant overlap of the pho- 
toionization and collisional ionization dominated regimes 
as was found in Section III A 3. 

Figures 6 and 7 lead to the following explanation of 
the sequence of events. During the laser pulse the cluster 
charges up and some outer-ionization takes place. After 
the laser pulse collisional ionization causes charge migra- 
tion from the outer shells to the core as electrons lose 
their energy by collisionally ionizing targets and falling 
deeper into the cluster's potential. The outer shells then 
explode faster than the inner shells which contain these 
cooled electrons. 



6. Kinetic energy 

The final kinetic energy distribution of the ions pro- 
vides information about how the cluster has evolved af- 
ter the laser pulse and is measurable in experiments. In 
Fig. 8 we plot the ion kinetic energy distribution (red 
solid line) for Ar2057, Ar56i, Ari47 and Ars5 clusters at 
the end of the simulations. To gain further insight into 
how the cluster disintegrates, we also plot the kinetic 
energy distributions for each charge state individually: 
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FIG. 7. (color online) Net charge per atom versus shell index 
for (a) Ar2057, (b) Ar 5 6i, (c) Ari47 and (d) Arss clusters after 
the laser pulse. The solid (blue) line is the average over all 
shells. 



(green) long dashes for Ar 1+ , (blue) medium dashes for 
Ar 2+ , (magenta) short dashes for Ar 3+ and (black) dash- 
dots for Ar 4+ . 

The largest clusters produce the most energetic ions. 
For all cluster sizes the Ar 2+ is the most energetic and 
dominates the high energy end of the spectrum. For each 
cluster size, there is an obvious "knee" in the total spec- 
trum that clearly follows the Ar 2+ spectrum. This is 
because the Ar 2+ is the most populous ion in all but the 
inner most shells. Note that the shell location of the 
ions was determined by detailed positional analysis (not 
shown here) and is not evident from Fig. 8. The Ar 2+ 
ions on the high energy side of the knee are from the 
exploding outermost shell and acquire a large amount of 
kinetic energy because there is a relatively low amount 
of electron screening. The Ar 2+ on the low energy side 
of the knee are from the other shells. 

Ar + appear as the least energetic ions. This is due to 
most of the Ar 4+ are produced in the core where collisions 
are more frequent. Since the core is more shielded than 
the outer shells, these Ar 4+ are effectively screened and 
thus gain less kinetic energy. 

The Ar 3+ populations decrease smoothly as a func- 
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FIG. 8. (color online) Ion kinetic energy distributions of (a) 
Ar2057, (b) Arsei, (c) Ari47 and (d) Arss clusters at the end 
of the simulations for each charge state. The solid (red) lines 
is the total ion spectrum. 



tion of kinetic energy for the larger clusters, and dom- 
inate the total spectrum at low energies (below around 
0.2 Hartree). There is a knee in the Ar 3+ spectrum for 
the smaller clusters, at 0.2 and 0.1 Hartree for Ari 4 7 and 
Arss, respectively, though the knee is somewhat washed 
out for the larger clusters. The ions with kinetic energy 
preceding the knee are in the core and expand slowly due 
to electron screening. The Ar 3+ dominates at low energy 
because it is the most populous ion in the core. 

The Ar 1+ ions are distributed almost evenly over the 
range of energy for the larger clusters with only a small 
peak at the high-energy end. In the smaller clusters, this 
peak is more pronounced. In all cases, this peak arises 
from ions on the outer most shell of the clusters, with 
energies about half of the energy of the most energetic 
Ar 2+ , indicating that the highest energy Ar 2+ and Ar 1+ 
are not screened. These high kinetic energy Ar 1+ (those 
at the peak) and Ar 2+ (those after the knee) primar- 
ily reside on the outermost shell and were ionized in the 
laser pulse solely by photoionization. They were thus not 
exposed to nor screened by low energy electrons. This ex- 



plains why there is no similar phenomenon of a large knee 
or peak for Ar 3+ as they must be created by collisional 
processes. 

In Fig. 8c) a significant decrease in the Ar 1+ spectrum 
is evident around 1 Hartree. This is consistent with the 
experimental observations in Ref. [ ] for similarly sized 
clusters, which stated that Ar 1+ fragments were mea- 
sured only up to 30 eV (1.1 Hartree). The largest contri- 
bution of high energy ions will come from the intensity 
peak of the laser pulse and from the largest clusters in the 
log-normal size distribution. The Ar^ are thus the most 
likely origin of these fragments and our model agrees well 
with the experimental ion kinetic energy observed. 

B. Electrons 

In the XUV, the laser interacts with the cluster 
through photoionization. The electrons are ejected from 
their parent ion or atom largely parallel to the laser polar- 
ization, but its subsequent motion is almost completely 
independent of the laser field due to the very small quiver 
energy at short wavelengths, even at these high intensi- 
ties. Thus electron motion is determined entirely by the 
Coulomb fields of the other charged particles in the clus- 
ter. After photonionization, electrons that do not have 
enough energy to escape the cluster space charge will 
disperse its energy throughout the cluster via collisions. 
The following sections provide details of how the energy 
is transferred, which requires a microscopic model of the 
relevant process to provide an accurate analysis of elec- 
tron energy distribution. 



1. Kinetic energy distribution 

The electron kinetic energy distribution for clusters of 
different sizes at different times is shown in Fig. 9. These 
times are: near the peak of the laser pulse at 30 fs (red 
solid line), shortly after the laser pulse at 90 fs (green 
large-dashed line), after a significant decrease in the in- 
tensity of the high energy tail at times indicated in the 
legends (blue medium-dashed line), and at the end of 
the simulation when the electron distribution no longer 
changes, again at times indicated in the legends (magenta 
short-dashed lines). 

The plots show that distributions at 30 fs are the only 
ones that deviate significantly from a single temperature 
Maxwellian distribution since the tails are not linear on 
the log plot. By 90 fs, they have become Maxwellian 
as indicated by the linear tails and this is verified by a 
quantitative analysis. Though not shown, the distribu- 
tions of each cluster size remain the same for at least 
another 60 fs, before the clusters begin exploding. Thus 
the electrons thermalize quickly, indicating a high degree 
of collisions. Once the clusters start to explode at longer 
times, the distributions become less energetic, as evident 
from the graphs. 
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FIG. 9. (color online) Electron kinetic energy distribution 
for (a) Ar 2 o57, (b) Arsei, (c) Ar M 7 and (d) Ar 55 . The time 
at which each distribution was calculated is indicated in the 
legends. 



At 90 fs, there are proportionally more low energy elec- 
trons in larger clusters compared with smaller clusters, 
which have a proportionally larger number of high energy 
electrons; the differences are about ~10 % between the 
Ar55 and Ar2057, so may be difficult to see directly from 
the graph. This indicates that in larger clusters, fast 
electrons generally lose their kinetic energy more read- 
ily than in small clusters. This is consistent with the 
finding of Section III A 2 that the larger clusters disperse 
the laser energy more quickly than small clusters. As 
there are more collisions in larger clusters, there will be 
faster thermalization and more collisional ionization, the 
latter of which actually removes energy from the ionized 
electron population. 

The aforementioned changes as the clusters evolve. At 
^150 fs (not shown) the distributions for all sizes are very 
similar, and at later times smaller clusters will have more 
lower energy electrons due to more rapid cluster explo- 
sion. Larger clusters have a larger quasi-neutral core and 
thus disintegrate more slowly; this allows for the preser- 
vation of more energetic electrons thus the distributions 
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for the larger clusters take a longer time to become less 
energetic. Note that the (blue) medium dashed line plots 
look very similar for each of the cluster sizes, but they 
were all calculated at different times, the latest time be- 
ing for the largest cluster. 



2. Velocity distribution 

The laser polarization direction sets the axis along 
which most photoionization occurs. Anisotropic pho- 
toelectron emission was observed in low intensity syn- 
chrotron experiments, though emission from clusters was 
less anisotropic than from atoms [3l ] . Only collisional 
processes can destroy the anisotropy inherent to pho- 
toionization. Since we have seen that clusters in intense 
laser pulses are highly collisional systems, which become 
more collisional as the size of the cluster is increased, we 
might expect that the effect of collisional processes would 
cause reduced anisotropy of the electron emission. 

Figure 10 shows the electron velocity distribution at 
different times for Ar2057, Ars6i, Ari47 and Arss clusters. 
The times chosen are the same as in the previous section, 
and are also listed in the legends of the figures. The dis- 
tributions parallel (lines) and perpendicular (symbols) to 
the laser polarization are shown for each cluster size and 
for each time. Though not shown explicitly, each distri- 
bution is symmetric about its respective axis. Data for 
both perpendicular directions are also included, though 
they are plotted with similar symbols of the same color 
for a given time. The degree of anisotropy is evident by 
the difference between the line and the points. 

At 30 fs, near the peak of the laser pulse, the low ve- 
locity electrons are increasingly isotropic as cluster size 
increases. This indicates that larger clusters thermalize 
more rapidly than the smaller clusters, and are thus more 
collisional. 

There is increasing isotropy at high velocity as cluster 
size increases at all times. The high velocity electrons are 
those which have escaped the cluster and are no longer 
subject to collisions. In the small clusters, they thus form 
a peak centered at 1.27 Hartree, the 3p photoelectron 
velocity. 

It can also be seen that in Ari47 clusters there is a 
larger decrease in intensity along the polarization axis 
from 30 fs to the final time at a velocity just below the 
peak, near 1 Hartree, than in Ar^. This is due to the 
Ai'147 clusters having a larger space charge and being 
more collisional. This results in a decrease of the veloc- 
ity of the later photoelectrons. The same is true in the 
larger clusters which are even more collisional and also 
have a larger space charge. Thus the Ar2057 clusters have 
the largest difference as seen in Fig. 10a. This trend ex- 
plains the lack of a clear 3p peak in the Ars6i and Ar2057 
clusters. 

The trend towards isotropy for larger clusters continues 
to increase for the highest velocities shown. The Ar 2 os7 
clusters show the same amount of anisotropy at 30 fs 
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FIG. 10. (color online) Electron velocity distribution parallel 
(lines) and perpendicular (symbols) to the laser polarization 
for clusters of size (a) Ar2057, (b) Ars6i, (c) Ari47 and (d) 
Arss. The time at which each distribution was calculated is 
indicated in the key. Note that both perpendicular directions 
are considered, and these are both plotted for each time with 
similar symbols of the same color, as indicated in the key. 
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as at the end of the simulation, indicating that many of 
the early photoelectrons have undergone some collisions 
before exiting the cluster. In smaller clusters, even more 
electrons leave the cluster unabated. In the Ar56i, cluster 
this causes a larger anisotropy gap between the 30 fs and 
the end of the simulation than in Ar2057- In the small- 
est clusters, this leads to the formation of the distinct 
photoelectron peaks. 

The anisotropy at early times observed in the small 
clusters for low velocity electrons has similar origin. 
There are a significant proportion of electrons which exit 
the cluster unabated. This leaves only a few remaining 
electrons to collisionally ionize or excite neutrals or ions, 
which losing some of their kinetic energy and thus veloc- 
ity. In larger clusters there are many more targets and a 
larger space charge at 30 fs which increases the number 
of collisions and removes the anisotropy at low velocity. 

The electron velocity distribution shows that larger 
clusters are more collisional and become so earlier. This 
leads to an increase in the isotropy of the emitted elec- 
trons. 



3. Connection with experiment 

Finally, we compare our simulations with experimental 
data of Ref. [2] by calculating the electron energy spec- 
trum over the intensity profile of the pulse. The spatial 
profile of the pulse was assumed to be Gaussian with a 
focus of 50 fim [2], and the cluster jet size was taken to be 
100 fim, the same size as the nozzle. The cluster size dis- 
tribution was assumed to be log-normal with (N) = 80 
and AN = 80, which was estimated via simulations with 
Ar§o as well as the two nearest closed-shell icosahcdral 
clusters, ie. Arss and Ari47 

We calculated the electron energy spectrum by inte- 
grating over cluster size and laser intensity profile, for 
intensities from 5 x 10 13 W/cm 2 to 1 x 10 12 W/cm 2 . 
Figure 11 shows our result (blue dashed line) compared 
with the experimental result (red solid line) for a peak 
intensity of 5 x 10 13 W/cm 2 [2]. We found reasonable 
agreement, given our lack of knowledge of the laser pro- 
file and the precise experimental setup. 

In Ref. [23] the role of augmented collisional ionization 
was shown to be necessary in order to obtain the highest 
charge states seen in Ref. [2], Ar 4+ . Figure 11 shows 
that the model is also capable of explaining the dominant 
features in the electron spectrum seen in the experiment. 
The success of this model is due to the more accurate 
treatment ionization processes. 



IV. SUMMARY 

We have shown the effect of cluster size on many as- 
pects of intense XUV-cluster interaction. Our model is 
verified by reasonable agreement with experimental ob- 
servations in Ref. [ ] . This includes reproducing the elec- 
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FIG. 11. (color online) The electron kinetic energy distribu- 
tion integrated over the spacial profile of a Gaussian pulse 
with a peak intensity of 5 x 10 13 W/cm 2 a focus of 50 /im and 
integrated over a log normal distribution of cluster sizes with 
(TV) = 80 and AN = 80 is shown as the (blue) dashed line. 
The experimental data (red solid line) is taken from Ref. [2], 



tron emission spectrum, obtaining the highest observed 
charge state, and obtaining a close match of the maxi- 
mum kinetic energy of the Ar 1+ species for Ari47 clusters. 

We find that, for all measures and aspects of the clus- 
ter, collisional processes contribute significantly. Further, 
the size of the cluster increases the importance of the col- 
lisional processes. Larger clusters proceed from predom- 
inantly photoionization driven to collisional ionization 
driven more rapidly than smaller clusters. This causes a 
significant modification in the total cluster photoioniza- 
tion cross-section versus what would be expected from a 
gas, since less neutrals are photoionized when collisional 
processes are included. This results in a decrease in the 
deposition of energy by the laser. Even for a small cluster 
with 55 atoms, neglecting collisional processes will over- 
estimate the amount of energy absorbed from the laser 
by over 30%. We term this process collisionally reduced 
photoabsorbtion (CRP). 

The charge states were shown to be in greater abun- 
dance, proportionally, in the larger clusters. The highest 
charge states appear earlier in larger clusters and their 
appearance happens well into the disintegration of the 
cluster. At the photon fluence considered, almost no 
charge states above Ar 3+ were observed during the laser 
pulse. Those that were were created almost exclusively 
by a two step collisional ionization process, wherein an 
ion is first collisionally excited, then ionized from the 
excited state, termed augmented collisional ionization 
(ACI). ACI was shown to be the dominant collisional 
channel in general. At any given time, around 20% of 
neutrals and ions are in an excited state. 

An examination of the charge of the cluster shells as 
a function of time show evidence of charge migration 
from the outer to inner shells. Immediately following 
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the laser pulse, the shells are almost uniformly charged, 
with the notable exception of the more highly charged 
outer shell in the larger clusters. However, as the system 
further evolves, collisions cause an increased charging of 
the outer shell of all cluster sizes, resulting the explosion 
of the outer shell and slow expansion of the inner shells, 
a process observed experimentally in mixed clusters [1]. 

An analysis of the ion kinetic energy distribution of 
each charge state species found that the high energy tail 
was almost entirely due to Ar 2+ . Most of the highest 
charge states were found to have very little kinetic en- 
ergy. This because they are created in core of the cluster 
where they are shielded by electrons during the cluster 
disintegration. 

The electron kinetic energy distribution was found to 
be close Maxwellian, except for very early times, indicat- 
ing rapid electron thermalization. The electron velocity 
distribution was found to be largely isotropic for large 
clusters, but highly anisotropic for small clusters. This 
is due to the high velocity electrons in the small clusters 
originating from photoionization events. The isotropy in 
the large clusters is evidence of high velocity electrons 



undergoing multiple collisions before escaping the clus- 
ter. 

Our findings may have implications for the direct 
single-shot imaging of large molecules with high inten- 
sity X-rays [? ]. These systems will be highly collisional, 
which may, for example, lead to a rapid change in the 
photoabsorbtion cross-section of the molecule via CRP. 
Models of this interaction will need to account for this 
effect. 
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